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Abstract 

We have studied odd-even nuclear mass staggering with the Skyrme-Hartree-Fock-Bogoliubov 
theory by employing isoscalar and isovector contact pairing interactions. By reproducing the 
empirical odd-even mass differences of the Sn isotopic chain, the strengths of pairing interactions 
are determined. The optimal strengths adjusted in this work can give better description of odd- 
even mass differences than that fitted by reproducing the experimental neutron pairing gap of 
120 Sn. 

PACS numbers: 21.60.Jz, 21.10.Dr 
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One of interesting phenomena in nuclei is the odd-even staggering (OES) of binding ener¬ 
gies. It is believed that OES is attributed to the pairing correlation, which plays an important 
role in nuclear structure [1, 2|. Numerous microscopic calculations, such as Hartree-Fock 


+ BOS (HF + BOS) or Hartree-Fock-Bogoliubov (HFB) theories, 
investigate the relationship between pairing interaction and OES 


rave been performed to 


13]. The interaction 


strength, Vo, of the pairing interaction is a crucial parameter to understand the nuclear 
properties of short-range correlations. 

Conventionally, the value of pairing strength is obtained by adjusting the average HFB 
pairing gap of even-even nuclei to fit the experimental odd-even mass differences of the 
neighboring nuclei. However, the OESs calculated by theoretical masses with the strength 
determined by the above pairing gap are in some cases substantial 


even mass differences of the experimental binding energies 
strengths adjusted to the average pairing gap in 120 Sn 14 


Hal- 


ly different from the odd- 
8 ]. The standard pairing 
were found to be too small to 
Hence, it is useful to refit the 


make a global comparison to the experimental OES data 
pairing strengths with global calculations to examine the whole system of the OES data. 

In the literatures, there are several measures of the empirical OES, such as three-point, 
four-point, and five-point formulae QS, 0 . Here we use the three-point formula defined 


as follows, 


A (3) = 


t^a+i 


[B(N + 1 ,Z)~ 2 B(N, Z ) + B(N - 1, Z)\, 


( 1 ) 


where B(N, Z) is the binding energy of the (A, Z) nucleus and 7 ta = (“l )" 4 is the number 
parity with A = A + Z. This second-order difference of binding energies is centered at an 
odd nucleus, i.e., odd-A nucleus for neutron OES. The A^ formula can reduce the mean- 
field contributions to the gap energy fi Q. For evenW nnelei. the OES is nrore sensitive to 
single-particle energies [3|, which is not discussed in this work. 

Our investigations are based on the self-consistent Hartree-Fock-Bogoliubov (HFB) cal¬ 
culations with Skyrme energy functionals in the particle-hole channel. We adopt the most 
commonly used Skyrme parameter sets, SLy4 15], SkP |i0, SkM* [i0, UNEDFO 10 and 
UNEDF1 []j|. In the particle-particle channel, we employ both isoscalar and isovector 
density-dependent delta pairing interactions. The isoscalar delta interaction is of the form, 


V( rr,r 2 ) = V 0 [l- ?/(—) 7 ]h(>i - r 2 ), 

Po 


( 2 ) 


where Vq is the pairing strength, rj and 7 are parameters, and p is the total density, while po 
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is the saturation density, which equals to 0.16 fm -3 . According to the choice of 77 , one can 
obtain different types of pairing, usually called volume, mixed, surface pairings. The volume 
interaction corresponds to 77 = 0 , which means that there is no explicit density dependence. 
It mainly acts inside the nuclear volume, while the surface pairing (77 = 1) is sensitive to 
the nuclear surface, and the mixed pairing (77 = 0.5) is a mix in these two pairings. In our 
calculations we choose 7 = 1 . 

The isospin-dependent pairing interactions have been proposed to reproduce better pair¬ 
ing gap in nuclei. One kind of isovector pairing, denoted by MSH pairing 
follows, 


6 ], is written as 


C“f r H (r,.r 2 ) = Vb[l - (1 - P)v,(-T- 

P 0 

-f3r}n(— )“ n ]5(ri - r 2 ), 
Po 


(3) 


where p = p n + p p , f3 = ( p n — p p )/p. The parameters are adjusted in the HFB framework 
to reproduce the pairing gaps in symmetric matter and neutron matter. Here, we adopt the 
best parametrization with p s = 0.598, a s = 0.551, p n = 0.947 and a n = 0.554 jf|. 


Another different isovector pairing, denoted by YS pairing 


9 ], is parameterized like this: 


^LvOr, r 2 ) = V 0 [1 - (770 + 77 ir 3 p){—) 

Po 

-V 2 W —) 2 ]<Kri - r 2 ), 

Po 


(4) 


0 . 




where 770 = 0.5, 771 = 0.2, 772 = 2.5, and 73 = — 1 for protons and 1 for neutrons 

We have carried out the HFB calculations with the latest version of HFBTHO |20j. The 
HFB solver HFBTHO has been developed by implementing the modified Broyden method 
and shared memory parallelism to accelerate the calculation speed. 

The even-even nuclei were first calculated in the HFB framework. We used the orbital 
space of 20 major harmonic oscillator shells, which is enough for the density functional 


calculations. For the pairing interactions, we adopted the cutoff energy of 60 
As for odd-A nuclei, we employed the equal filling approximation (EFA) 


MeV. 


21]. Starting 


from the HFB solution of neighboring even-even nuclei, we select quasi-particle orbitals 
for the blocking of the odd nucleon. The one-quasiparticle configurations are determined 
within the blocking energy window A ]qpwin , which is smaller than 8 MeV for light nuclei and 
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bigger than 2 MeV for heavy nuclei. In present work, we took £d qP]W i n = 2h/\fA. Finally, 
we performed unconstrained self-consistent calculations for all candidate configurations and 
took the minimum energy as the binding energy of the odd-A nucleus. 


TABLE I. The optimal pairing strengths Vo adjusted to give the best fit to the neutron odd-even 
staggering of the Sn isotopic chain. The word ’standard’ means the pairing strengths fitted by the 
average pairing gap of 120 Sn. The unit of the pairing strengths is MeV-fm 3 . 



SLy4 

SkP 

SkM* 

UNEDF0 

UNDEF1 

Standard 

283.3 

213.3 

233.9 

- 

- 

Mixed 

310 

240 

270 

225 

245 

MSH 

400 

325 

360 

300 

325 

YS 

325 

260 

290 

245 

260 


As mentioned before, the pairing strengths should be refitted to examine the whole system 
of the OES data. In principle, global calculations are needed. However, such calculations for 
different Skyrrne forces and pairing interactions are very time consuming. We only adjust the 
pairing strengths to give the best fit to the OES of the semimagic Sn isotopes with neutron 
number ranging from 49 to 85. In fact, Sn isotopes are known to be excellent laboratories for 
comparison with mean field calculations and OES effects 10]. In present work, the proton 
pairing strengths are identical to the neutron ones. 

Table I lists the optimal pairing strengths for the five different Skyrrne forces with different 
pairing interactions. The results for the root mean square deviation of our calculations are 
shown in Figure 1. The root mean square deviation a of OES between HFB calculations 
and experimental data is defined as 


cx(OES) = 


Ni 


^|A< 3 >(HFB)-A j (3 >(Exp)|VJV i 

\ <=1 


(5) 


where N t is the number of data points. The experimental masses are taken from Ref. 22], 


In Figure 1, we compare the results for the three commonly used Skyrrne parameter sets, 
SLy4, SkP and SkM*. It seems that the optimal pairing strength to obtain the best fit to 
binding energies are different from that for OES. The minima of the fitting curve of binding 
energies for SLy4 and SkM* forces with mixed pairing are slightly bigger than that of OES, 
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FIG. 1. (Color online) Root mean square deviation o between HFB calculations and experimental 
data. The upper, middle and lower panels correspond to HFB calculations with Mixed pairing, 
MSH pairing and YS pairing, respectively. Filled symbols stand for the mean square deviation of 
binding energies of the Sn isotopic chain, while open ones for the deviation of OES. The gray stars 
in the upper panel represent the results calculated with the strength fitted by the average pairing 
gap of 120 Sn. See text for details. 
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FIG. 2. (Color online) Calculated and experimental values of Ao ; for the semimagic Sn and Pb 
isotopic chains. The SLy4 interaction is adopted with mixed pairing in the HFB framework for 
different pairing strengths. The dark cyan down triangles shows the results with the standard 
pairing strength listed in Table 1 and the blue up triangles stand for the results with the pairing 
strength fitted by Bertsch et. al. [8(] 


while for the SkP force, it shows an opposite result. The small differences of the minima 
for binding energies and OES come from the definition of the root mean square deviation 
cr. From Eqs.(l) and (5), we can see that the deviation a of OES contains the correlations 
between neighborhood nuclei, which is not included in the deviation of binding energies. 

For the SkP force, the different results from SLy4 and SkM* forces may come from its 
parametrization procedure. The SkP parameters are optimized together with pairing force, 


which is similar to the recently published parameter sets UNEDF0 and UNEDF1 


5 , 118 , 


19]. 


In fact, we also adjusted the pairing strengths for UNEDF0 and UNEDF1 forces with mixed 
pairing interaction. The results show that as the pairing strength increases, the root mean 
square deviation of binding energies increases sharply. Take UNEDFO+Mixed as an example, 
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FIG. 3. (Color online) Calculated and experimental values of Ao ; for neutrons in the Sn and 
Pb semimagic isotopic chains. SLy4 is adopted together with the mixed pairing (Eq.(2)) and two 
isospin dependent pairing (MSH in Eq.(3) and YS in Eq.(4)) interactions in the HFB model. 


the optimal pairing strength for OES is at 225 MeV • fm 3 , where the mean square deviation 
cr of binding energies is 3.18 MeV. We remind the reader that the fitting of pairing strengths 
should be done with caution for some Skyrrne forces. 

Figure 2 shows the OESs of two semimagic isotopic chains Sn and Pb with three pairing 
strengths fitted by different ways. The overall trend is reproduced for all three treatments of 
pairing. For Sn isptopes, the flatness up to the quenched gap at N = 83 is well consistent to 
experimental data, and for Pb isotopes, the downsloping trend up to gap at N = 125 is also 
reproduced. The pairing strength fitted by overall systematics by Bertsch et. al.js| gives 
higher average OES in both spherical chains, while the strength adjusted by reproducing the 
pairing gap of 120 Sn gives too small values in the Sn isotopic chain. Our results are between 
these two values. 

We have also compared the results calculated by the optimal strengths with three different 


7 









pairing interactions, shown in Figure 3. The MSH and YS pairing interactions are essentially 
the mixed-type contact pairing interactions with isospin dependence. The overall trends for 
the three different pairing interactions are similar. For nuclei with neutron excess, the OESs 
with MSH and YS pairing interactions are smaller than with mixed pairing interaction. The 
isospin dependent pairing interactions flatten the odd-even mass differences as a function of 
neutron number, which is consistent with the results in Ref. [7]. 

In summary, we have investigated the neutron OESs of Sn and Pb isotopes using self- 
consistent Skyrme-Hartree-Fock-Bogoliubov theory with SLy4, SkP, SkM*, UNEDFO and 
UNEDF1 forces together with mixed pairing and two different isospin dependent pairing 
interactions. The pairing strengths are adjusted by reproducing the empirical OESs of Sn 
isotopes. The pairing strengths necessary to obtain the best fit to the binding energies are 
different from that for OES. We reproduced the flatness of the OES due to the isospin effects 
and compared to the results with isoscalar pairing interactions. 
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